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Abstract 


This paper examines the impact of geoengineering via stratospheric sulfate 
aerosol on the quasi-biennial oscillation (QBO) using the NASA Goddard Earth 
Observing System (GEOS-5) Chemistry Climate Model. We performed four 30- 
year simulations with a continuous injection of sulfur dioxide on the equator at 
0° longitude. The four simulations differ by the amount of sulfur dioxide injected 
(5Tg/year and 2.5 Tg/year) and the altitude of the injection (16km-25km and 
22km-25km). We find that such an injection dramatically alters the quasi- 
biennial oscillation, prolonging the phase of easterly shear with respect to the 
control simulation. This is caused by the increased aerosol heating, and 
associated warming in the tropical lower stratosphere and higher residual 
vertical velocity. In the case of maximum perturbation, i.e. highest stratospheric 
aerosol burden, the lower tropical stratosphere is locked into a permanent 
westerly QBO phase. 
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Introduction 


The quasi-biennial oscillation (QBO) is an approximately 28-month period 
oscillation of zonally symmetric easterly and westerly winds in the tropical 
stratosphere. The QBO is caused by vertically propagating waves, such as 
equatorial Kelvin and Rossby-gravity waves, that deposit momentum in the 
stratosphere [e.g. Baldwin et al, 2001], While the QBO is confined to the tropics, 
its phase affects the stratospheric transport to the extratropics and the strength 
of the polar vortex [Holton and Tan, 1980], altering transport from the tropics to 
mid- and high latitudes of stratospheric trace gases and aerosols [ Trepte and 
Hitchman, 1992], Additionally, studies have shown that the QBO can impact the 
tropospheric winds ( Garfinkel and Hartmann, 2011) and precipitation (Jihoon et 
al. 2013). 

The vertical descent of the QBO wind shear is linked to the mean tropical 
upwelling of the Brewer-Dobson circulation (BDC) [e.g. Watanabe and Kawatani, 
2012], Kawatani and Hamilton (2013) identified in radiosonde observations for 
the 1953-2012 period a long-term trend of weakening QBO amplitude, which 
they attributed to the BDC strengthening due to increasing greenhouse gases. 
Geoengineering is a deliberate modification of the Earth system in order to 
counteract global warming due to increasing greenhouse gases. Some proposed 
geoengineering methods address the causes of the warming by reducing the 
amount of solar radiation reaching the Earth surface (solar radiation 
management or SRM), for instance by continuously injecting sulfate aerosol into 
the stratosphere. This method aims to reproduce the global surface cooling 
observed after major volcanic eruptions. 
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In addition to its effect on tropospheric temperatures ( Kravitz et al., 2013), 
precipitation ( Haywood et al., 2013; Tilmes et al., 2013) and stratospheric ozone 
( Tilmes et al., 2009), such stratospheric injection of aerosol would also lead to 
perturbations of stratospheric dynamics. An increase in stratospheric aerosol 
loading would warm the lower stratosphere, mainly via absorption of longwave 
radiation. Such warming would lead to a strengthening of the tropical upwelling, 
as showed by Aquila et al. (2012) in the case of a Mt. Pinatubo-like eruption, 
which could interfere with the periodicity of the QBO. After the Mt. Pinatubo 
eruption, observations showed a warming of the lower stratosphere of about 3K, 
and a delay in the downward propagation of the easterly shear [ Labitzke , 1993], 
Here, we present a set of four model experiments that simulate the 
geoengineering stratospheric injection of sulfur dioxide, varying the burden and 
altitude of the SO 2 injection, in order to investigate the impact of stratospheric 
aerosol geoengineering on the QBO. 

Model simulations 

The Goddard Earth Observing System Chemistry Climate Model (GEOSCCM) 
couples the GEOS-5 general circulation model [Rienecker et al., 2011], the 
Georgia Institute of Technology-Goddard Global Ozone Chemistry Aerosol 
Radiation and Transport (GOCART) module [Colarco et al., 2010] and the 
StratChem stratospheric chemistry module [Pawson et al, 2008], GOCART is a 
bulk aerosol model which include a parameterization of the chemical production 
of SO 4 aerosol from oxidation of dimethyl sulfide (DMS) by OH during day and 
NO 3 during night, and from oxidation of SO 2 by OH in the gas phase and by 
H 2 O 2 in the aqueous phase. Evaluation of a similar version of GEOSCCM with 
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respect to stratospheric aerosol by Aquila et al. [2012; 2013] shows good 
agreement with observations of aerosol distributions and ozone and NO 2 
depletion after the eruption of Mount Pinatubo. Compared to Aquila etal. (2013), 
this GEOSCCM version includes both a mechanism to generate the QBO using a 
gravity wave drag parameterization [Molod et al., 2012], and a coupling between 
aerosol and heterogeneous chemistry via the aerosol surface area density. The 
aerosol surface area density is calculated from the dry sulfate mass assuming 
that the aerosol particles are lognormally distributed with modal radius 0.35 pm. 
This same size distribution, hydrated accordingly to the relative humidity, is 
used to calculate the optical properties of the stratospheric sulfate aerosol and 
its settling velocity. 

The GEOSCCM resolution is 2.0° latitude by 2.5° longitude, with 72 vertical 
hybrid levels from surface to 0.01 hPa. The model is prescribed with sea surface 
temperatures and sea ice concentrations calculated with the Community Earth 
System Model (CESM) [Gent et al., 2011] using emission inventories valid for the 
Representative Concentration Pathway (RCP) 4.5 [Taylor et al., 2012], 

We performed four 30-year long experiments from 2020 to 2049 in which we 
prescribed a continuous injection of SO 2 in the stratosphere on the equator at 0 
longitude. The four perturbed simulations differ from each other with regard to 
the burden (5 Tg/year in two of the experiments and 2.5 Tg/year in the other 
two) and altitude (16km-25km and 22km-25km) of the geoengineering injection. 
The four experiments (G 5 16 " 251 ™, G 2 . 5 16 25km , Gs 22 25km and G 2 . 5 22_25km ) are 
summarized in Table 1. Gs 16 25km corresponds to the experiment G4 of the 
Geoengineering Model Intercomparison Project (GeoMlP; Kravitz et al., 2011). 
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Additionally, we performed a control simulation without stratospheric SO 2 
injection. 

We concentrate in our discussion on the last 20 years of simulations, when the 
sources and sinks of geoengineering aerosol are in equilibrium. During this time 
span, the atmospheric burden of the geoengineering sulfate aerosol is equal to 
4.7 Tg-S in G 5 22 ' 25km , 3.1 Tg-S in Gs 16 ' 251 ™, 2.1 Tg-S in G 2 . 5 22 ~ 25km and 1.5 Tg-S in 
G 2 . 5 16 - 25km (Table 1). The mixing between the tropics and extra-tropics is weaker 
at the altitudes where the aerosol is injected in Gs 22 25km and G 2 . 5 22_25km than at 
lower altitudes (16-20 km), and this leads to longer stratospheric residence 
times and higher aerosol burdens in these experiments (Table 1). For 
comparison, Baran and Foot (1994) measured a maximum stratospheric aerosol 
burden of about 7 Tg-S after the eruption of Mt. Pinatubo, which decreased to 4 
Tg-S about 18 months after the eruption. The vertical distribution of the 
geoengineering aerosol in the four simulations is shown in Fig. 1 of the 
supplementary material. The total burden in our simulations are similar to that 
in the geoengineering experiments of Heckendorn et al. (2009), but the aerosols 
are more confined to the tropics in our simulations. This is likely due to 
differences in the transport (Heckendorn et al. (2009) use a two- 
dimensional model), but could also be due to differences in the injection height 
and microphysical processes. 

The stratospheric aerosol optical thickness (AOT) can be roughly converted to 
radiative forcing by multiplying by 25 W/m 2 ( Hansen et al, 2005). Using this 
conversion factor, the aerosol radiative forcing in our simulations ranges from - 
lWnr 2 in G 2 . 5 16 " 25km to -3 Wnr 2 in G 5 22 " 251 ™ 1 (see Fig. 2 supplementary material for 
the simulated AOT). Even in a fairly optimistic scenario such as the 


©2014 American Geophysical Union. All rights reserved. 





Representative Concentration Pathways (RCP) 4.5 [ Clarke et al., 2007; Smith et 
al., 2006; Wise et al, 2009], the radiative forcing by increasing greenhouse gases 
reaches 3Wnr 2 around 2030. Hence, the aerosol perturbations introduced in our 
experiments are within the possible range that might be needed to offset 
warming from increasing greenhouse gases. 

Note that in these experiments, the source function for gravity waves is held 
fixed in all experiments and is not coupled to convection, and the sea surface 
temperatures (SSTs) do not vary among the experiments. It is likely that the SSTs 
and the convection will change were geo-engineering implemented. Hence, the 
source of wave driving for the QBO in these experiments does not vary 
realistically in response to geo-engineering, and future work is necessary to 
explore whether these limitations affect the model results presented below. 
Results 

Our results show that the stratospheric aerosol injection dramatically perturbs 
the QBO periodicity, prolonging the phase of easterly shear with respect to the 
control case. Figure 1 shows the vertical profiles of the tropical zonal wind, 
averaged between 2°S and 2°N, in the control simulation (upper panel) and in the 
four geoengineering experiments. The phase of easterly shear persists for longer 
and longer with increasing stratospheric burden of aerosol. Table 1 reports the 
mean value and standard deviation of the simulated QBO period: the period of 
the QBO increases with the burden of geoengineering aerosol, from about 25 
months in the control simulation to about 50 months in Gs 16 " 25km . The QBO 
completely disappears in Gs 22 25km (lowest panel), where the stratosphere is in a 
perpetual easterly shear or westerly phase. In the following discussion we will 
focus on G 5 22 25km . The response of the QBO to the increase of aerosol is fast, 
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within one QBO period in the G 5 experiments, and within two QBO periods in 
G 2 . 5 . Additionally, between 17 km and 30 km altitude the amplitude of the QBO 
decreases with increasing stratospheric aerosol burden (Fig. 3 of supplementary 
material). 

Figure 2 shows the mechanism that leads to the interruption of the QBO in G 5 22 " 
25km . The heating from the geoengineering aerosol induces a warming of the 
lower stratosphere up to about 27 km altitude (Fig. 2, upper left panel). Because 
the radiative damping timescales are longer in the lower stratosphere than in the 
upper stratosphere, the aerosols affect lower stratospheric temperature more 
strongly (Newman and Rosenfield, 1997). The warmer temperature centered on 
the equator induces an equatorial positive wind shear [J3u z = -RT /H), which 

results in persistent westerly winds between 20 and 30 km. The heating anomaly 
due to the aerosol also causes an increase of the residual vertical velocity w* 
(Fig. 2, upper right panel), which advects the aerosol upward from the initial 
injection altitude. The increase in w* extends well above the level at which the 
diabatic heating ends, consistent with Holloway and Neelin (2007) and Garfinkel 
etal. (2013). 

The typical processes that lead to the QBO's downward propagation weaken in 
the presence of the increase to w'due to the geoengineering aerosols. This effect 
is demonstrated by the lower panels of Fig. 2, which show the budget of the 
terms that force the QBO in the control and G5 22 " 251 ™ experiments, averaged 
between 2S and 2°N. When the QBO propagates downwards, momentum 
deposition from the parameterized gravity waves (GWD) drives the downward 
propagation, but it is opposed (but not fully) by westerly momentum from the 
vertical advection of the mean flow (-w*u z ) which advects the wind anomalies 
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upwards. In the presence of aerosols and enhanced w*, the easterly momentum 
from the parameterized gravity waves (GWD) and the westerly momentum from 
the vertical advection of the mean flow (-w*u z ) balance each other leading to a 
smaller time derivative of the zonal wind. 

Figure 3 shows the anomalies of temperature, zonal winds, and ozone and N 2 O 
concentrations in Gs 22 25km with respect to the control simulation. The circulation 
anomaly is also superimposed on Fig. 3. The temperature anomalies (Fig. 3, 
upper left) shows the largest warming centered at about 20 km and related to 
the absorption of longwave radiation by the aerosol. The cold anomaly above this 
region near 30km is related to the circulation anomaly caused by the aerosol. 
This induced circulation is also clearly reflected in the N 2 O concentrations (Fig. 
3, lower right). N 2 O is a well-suited tracer to study stratospheric transport, due 
to its long stratospheric lifetime and its distribution with higher concentrations 
at the surface in the tropics. The circulation changes induced by the aerosol 
increase the N 2 O concentrations in the middle stratosphere and enhance its 
transport to the extratropics. The ozone anomalies (Fig. 3, lower left) are due to a 
combination of heterogeneous chemistry on the aerosol particles and to the 
induced change in tropical dynamics ( Aquila et al. 2013), which advects air with 
different ozone concentrations. Similar changes are also found in the other 
simulations (see supplementary material); there are some differences in the 
spatial structure of the anomalies among the experiments, and future work is 
needed to examine these differences 
Conclusions 

Our simulations show that geoengineering injection of stratospheric aerosol can 
lead to dramatic changes in the QBO, prolonging the phase of easterly shear with 
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respect to the control simulation. For very large increases in stratospheric 
aerosol burden (4.7 Tg-S), the lower tropical stratosphere is locked into a 
permanent westerly QBO phase. Fig. 7 in the supplementary material suggests 
that there is roughly a quadratic relationship between the aerosol burden and 
the QBO period, calculated at 30 hPa. This modification of the QBO occurs 
because the increase in aerosol burden leads to a warming of the tropical lower- 
middle stratosphere, mainly via absorption of longwave radiation, and, hence, to 
stronger westerly winds resulting from the thermal wind relation. This warming 
also induces an increase in the residual vertical velocity w*, which broadens the 
band of westerly winds and lofts the aerosol higher in the stratosphere. Because 
the QBO can impact stratospheric and tropospheric ozone (Randel and Wu, 1996; 
Ziemke and Chandra, 1999; Oman etal, 2013), tropospheric winds ( Garfinkel and 
Hartmann, 2011), and precipitation (Jihoon et al. 2013), this geoengineering- 
forced permanent lower stratospheric QBO westerly phase could substantially 
alter surface climate. 
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Table 1: Summary of the model experiments performed for this study. The 


period of the QBO is calculated at 30 hPa. 
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Fig-1 : Vertical profiles of the zonal wind, zonally averaged between 2°S and 2°N 
in (from top to bottom) control run, G 2 . 5 16 25km , G 2 . 5 22_25km , Gs 16 251 ™, G5 22 25km . 
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Fig.2: Upper panels: Vertical profiles of the zonal mean (2S -2°N) anomaly of the 
temperature anomaly (upper left, in K) and residual vertical velocity (upper 
right, in mm/s). The anomalies are calculated between the two months of the 
control and Gs 22 25km when the u wind profiles are the most similar. The dashed 
blue line shows the aerosol heating rates in K/day in G 5 22 25km (note the scaling 
factor in the upper left figure). The dotted black line in the upper right panel 
shows iv’ in the control simulation. Lower panels: Vertical profiles of the forcing 
terms (in m s 1 day 1 ) that drive the QBO in the control (left) and G 5 22 25km (right) 
simulations during the months considered in the upper panels. Shown are the 
2°S-2°N zonal means of u wind in m/s multiplied by a factor 0.03 (dashed black), 
the time derivative of u (solid black), the advection of the mean flow (red), the 
parameterized gravity wave drag (green) and the divergence of the Eliassen- 
Palm flux (blue). 
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Fig-3 : Vertical profiles of the zonal mean anomalies of temperature, zonal wind u, 
ozone and N 2 O concentrations in G 5 22 25km with respect to the control simulation. 
The streamlines show the anomaly of the residual circulation. The anomalies are 
calculated over the whole 2030-2050 period. The thickness of the streamlines is 
proportional to the change of speed. 
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